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Mixtures of relativistic gases are analyzed within the framework of Boltzmann equation by using 
Grad's moment method. A relativistic mixture of r constituent is characterized by the moments 
of the distribution function: particle four-flows, energy-momentum tensors and third-order moment 
tensors. By using Eckart's decomposition and introducing 13r + 1 scalar fields - related with the 
four-velocity, temperature of the mixture, particle number densities, diffusion fluxes, non-equilibrium 
pressures, heat fluxes and pressure deviator tensors - Grad's distribution functions are obtained. 
Grad's distribution functions are used to determine the third-order tensors and their production 
terms for mixtures whose constituent's rest masses are not too disparate, so that it follows a system 
of 13r + 1 scalar field equations. By restricting to a binary mixture characterized by the six fields of 
partial particle number densities, four-velocity and temperature, the remainder 21 scalar equations 
are used to determine the constitutive equations for the non-equilibrium pressures, diffusion fluxes, 
pressure deviator tensors and heat fluxes. Hence the Navier-Stokes and generalized Fourier and 
Fick laws are obtained and the transport coefficients of bulk and shear viscosities, thermal conduc- 
tivity, diffusion, thermal-diffusion and diffusion-thermal are determined. Analytic expressions for 
the transport coefficients in the non-relativistic and ultra-relativistic limiting case are given. Fur- 
thermore, solutions of the relativistic field equations for the binary mixture are obtained in form of 
forced and free waves. In the low frequency limiting case the phase velocity and the attenuation 
coefficient are determined for forced waves. In the small wavenumber limiting case it is shown that 
there exist four longitudinal eigenmodes, two of them corresponding to propagating sound modes 
and two associated with non-propagating diffusive modes. 
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I. INTRODUCTION 



We may state that the beginning of the relativistic kinetic theory goes back to 1911 when Jiittner [1] derived 
an equilibrium distribution function for a relativistic gas, which in the non-relativistic limiting case becomes the 
Maxwellian distribution function. Jiittner has also succeeded to derive in 1928 [2] the relativistic Bose-Einstein 
and Fcrmi-Dirac distribution functions. The covariant formulation of the Boltzmann equation was proposed by 
Lichncrowicz and Marrot [3] in the forties of the last century and in the sixties the determination of the transport 
coefficients from the Boltzmann equation by using the Chapman-Enskog methodology was obtained by Israel [4] and 
Kelly [5] . 

Mixtures of relativistic gases are important in the field of astrophysics, in particular to problems associated with 
gases at high temperature in a stellar interior. Within the framework of Boltzmann equation these mixtures were 
studied by several authors and among others we quote the works [6-18]. In these works the usual methodology applied 
to determine the transport coefficients was the Chapman-Enskog method and general expressions for the transport 
coefficients were determined for generic differential cross sections. Only few explicit expressions for the transport 
coefficients were determined for very special cases of binary mixtures of Maxwellian particles [8] and hard spheres [13] 
where the constituents have the same rest masses or a Lorentz gas of hard spheres [14] . 

In this work we analyze mixtures of relativistic gases of Maxwellian particles by using Grad's moment method 
applied to the system of Boltzmann equations. In relativistic kinetic theory there is no unique differential cross 
section which tends in the non-relativistic limiting case to the one which is proportional to the relative velocity, 
known as the differential cross section of Maxwellian particles. At least three of them given in the works [4, 19, 20] 
have this property, and here we shall use the one defined in [20] . 

We are interested in a mixture that is characterized by r constituents and described by 13r + 1 basic fields of four- 
velocity, temperature of the mixture, particle number densities, diffusion fluxes, non-equilibrium pressures, heat fluxes 
and pressure deviator tensors. Grad's distribution functions are obtained from the definition of the basic fields which 
are related with the Eckart decomposition of the particle four-flows and of the energy-momentum tensors. From the 
knowledge of Grad's distribution functions the third-order moments of the distribution functions are calculated and 
the production terms that appear in the balance equations arc determined for the case of a mixture of Maxwellian 
particles for constituents whose rest masses are not too disparate. Once the 13?* + 1 field equations for the basic 
fields are established, the mixture is restricted to a binary mixture characterized by six fields of particle number 
densities, four-velocity and temperature of the mixture. The constitutive equations for the diffusion flux, pressure 
deviator tensors, non-equilibrium tensors and heat fluxes are determined by a method akin to the Maxwellian iteration 
procedure applied to the remaining 21 scalar field equations. The generalized laws of Fick, Fourier and Navier-Stokcs 
are obtained and the coefficients of diffusion, thermal-diffusion, diffusion-thermal, thermal conductivity and bulk and 
shear viscosities are determined. Explicit expressions for these coefficients in the non-rclativistic and ultra-relativistic 
limiting cases are given, as well as their graphics as functions of a parameter which represents the ratio of the rest 
energy of a particle and the thermal energy of the gas. From the system of field equations for the binary mixture 
we have analyzed the solutions corresponding to small perturbations from an equilibrium state related with the 
propagation of forced and free waves. For the acoustic solution in the low frequency limiting case we have determined 
the phase velocity and the attenuation coefficient. For the eigenmodes in the small wavenumber limiting case we have 
shown the occurrence of two propagating sound modes and two non-propagating diffusive modes. 

This work is structured as follows: in Section II the Boltzmann equations, the moments of the distribution functions 
and their balance equations are introduced. The Eckart decomposition is the subject of Section III and Grad distribu- 
tion functions are determined in Section IV. In Section V the constitutive equations for the third-order moment tensors 
and for the production terms are obtained and in Section VI the linearized 13r + 1 field equations are established. 
A binary mixture with six scalar fields of particle number densities, four-velocity and temperature of the mixture is 
analyzed in Section VII, where the generalized laws of Fick, Fourier and Navier-Stokes are obtained from a method 
akin to the Maxwellian iteration procedure. Furthermore, the transport coefficients of diffusion, thermal-diffusion, 
diffusion-thermal, thermal conductivity and bulk and shear viscosities associated with these laws are determined in 
the non-relativistic and ultra-relativistic limiting cases. In Section VIII solutions of the binary mixture field equations 
are analyzed in terms of sound propagation and eigenmodes. Finally, the main conclusions of this work are discussed 
in Section IX. 

Latin indexes running from a = 1 . . . r specify the constituents of the mixture, while Greek indexes running from 
a = 0, 1, 2, 3 denote the space-time coordinates. 
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II. BOLTZMANN AND TRANSFER EQUATIONS 

Let us consider a relativistic gas mixture of r constituents in a Minkowski space with metric tensor {rj a p) = 
diag(l, —1,-1,-1). The particles of constituent a = 1, . . . r have rest mass m a and are characterized by the space- 
time coordinates (x a ) = (ct,x) and momenta (p") = (p°,p Q ). The length of the momentum four- vector is m a c so 
that p° a = v/|p Q | 2 + m 2 a c 2 . 

An elastic collision between the particles of two constituent a and b whose momentum four- vectors before collision 
are denoted by p" and p" reads 

Pa+Pb^Pa+Pb, (2.1) 

where the primes denote the momentum four-vectors after collision. 

A state of the relativistic mixture of r constituents is characterized by the set of one-particle distribution functions 

/(x,p OJ t) = / OJ a = l,2,...,r (2.2) 

such that /(x, p a , t)d 3 xd 3 p a gives at time t, the number of particles of constituent a in the volume element d 3 x about 
x and with momenta in the range d 3 p a about p a . 

The one-particle distribution function of constituent a satisfies a Boltzmann equation, which in the absence of 
external fields reads (see e.g. [16, 17]) 

d 3 p b 



Pa dafa = T [ tf'Jh ~ fafb)F ba a ab , (2.3) 



where a ab and d£l denote the invariant differential clastic cross-section and the element of solid angle that characterizes 
a binary collision between the particles of constituent a with those of constituent 6, respectively. Moreover, F ba is the 
invariant flux defined by 



F ba = P y b] j (g - |) 2 - (g x |) 2 = yj (p >ba y - m W b c\ (2.4) 

and the following abbreviations were introduced 

/£ = /(x,p' a> *)> K = f(^Plt), / a = /(x,p a ,<), / 6 = /(x,p b ,t). (2.5) 

The general equation of transfer for the constituent a of the mixture is obtained through the multiplication of the 
Boltzmann equation (2.3) by an arbitrary function of the momentum four-vector ip a = ipip") and integration of the 
resulting equation over all values of d 3 p a /p a o, yielding: 



da f i>Ma— = E / Mf'Jl - fafb)F ba a ab dn 
J PaO ~ J 



d 3 p b d 3 p a 

PbO PaO 



/d^ pi) d^ p 
{lp' a - lpa)fafbFbaO-abdtt ° 
PbO PaO 



(2.6) 



where the last equality on the right-hand side of the above equation follows by using the symmetry properties of the 
collision term. 

If we sum (2.6) over all constituents we get a general equation of transfer for the mixture that reads 

3«£ / taPU*— = \ E A< + ^ ~ ^ - — — • ( 2 - 7 ) 

f^J PaO p ba PaO 

The right-hand side of (2.7) was also obtained by using the symmetry properties of the collision term and by inter- 
changing the dummy indexes a and b in the sums. 

The moments of the distribution function we are interested are the partial particle four-flow TV", the partial energy- 
momentum tensor T"^ and the partial third-order moment tensor T"^ 7 which are defined in terms of the distribution 
function through: 

K = c(p-J a ^, T? = cfptfjJ^, Tf^cfpyUfa^- (2.8) 



PaO J PaO J PaO 
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The particle four-flow of the mixture N a , the energy-momentum tensor T a/3 and the third-order moment T Q ^ 7 of the 
mixture are obtained by summing the partial quantities over all constituents, i.e. 

r r r 

N a = J2 N a> T al3 = T a' 3 , T afil = T a Pl - ( 2 - 9 ) 

a—1 a— 1 a— 1 

The balance equations for the partial particle four-flow, for the partial energy-momentum tensor and for the partial 
third-order moment tensor are obtained by choosing ip a = c, ip a = cp@ and ipa = cPaPa m (2-6), yielding 

d a N« = 0, d a Tf = ?l d a Tf* = ?l\ (2.10) 

respectively. Above, the production terms Pf and Pf 7 are given by 

pf = X> / (tf-Muf^v^i***, (2 . u) 

r-f J PbO PaO 

o—l 

P a 7 = E c / Vifl -P^2)faf b F ba a ab dQ^^. (2.12) 

r~f J PbO PaO 



6=1 



By summing (2.10) over all constituents of the mixture we get the balance equations for the particle four-flow, 
energy-momentum tensor and third-order moment tensor of the mixture that read 

r 

d a N a = 0, d a T al3 = 0, d a T a ^ = P M = E Pa 7 - (2-13) 

a=l 

These are the balance equations used in the relativistic theory of extended thermodynamics of a single fluid (see [21]). 

III. ECKART DECOMPOSITION 

For the decomposition of Eckart we introduce the four- velocity U a such that U a U a = c 2 and the projector 

A Q/3 = if 13 — \u a ljP, with A al3 U a = Q. (3.1) 

Hence, we write the partial particle four-flow as 

N% = n a U a + J", such that l%U a = 0. (3.2) 

Here n a is the particle number density and J" the diffusion flux of the constituent a in the mixture. They are defined 
in terms of the following projections of the partial particle four-flow 

n Q = ^N°U a , r a = A<$N? . (3.3) 

From the definitions (2.8a) and (3.2) we can represent the diffusion flux as 

J« = A£c/^/ ^. (3.4) 

J PaO 

Further by summing (3.2) over all constituents it follows 

r r 
a—1 a—1 

which gives the particle number density of the mixture as the sum of the partial particle number densities and shows 
- due to the constraint (3.5b) - that there exist only (r — 1) partial diffusion fluxes that are linearly independent for 
a mixture of r constituents. 

The decomposition of the partial energy-momentum tensor in the Eckart description is written as (see e.g. [22-24]) 
Tf = P<^ - (P Q + ^aW + (qf + h B jf) + Itf" K + h a J«) + ^U a U^ (3.6) 
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where we have introduced the partial quantities: p a hydrostatic pressure, w a non-equilibrium pressure, e a energy per 
particle, q" heat flux, h a = e a + p a /n a enthalpy per particle and pi"^ pressure deviator tensor, i.e., the traceless part 
of the pressure tensor. These fields are given in terms of the projections of the partial energy-momentum tensor by: 



f 



p a + K7 a = --A Q/3 T° 



q° + hj« = A«TfU p , p<^> = (a?A| - l^ aP ^ T2 S . 



(3.7) 
(3.8) 



Note that the symmetric and traceless part of a second order tensor Y a @ is defined in terms of the projector A a ^ 
through the relationship 



(a°a£ + A£A^ 



A"^A 70 



(3.9) 



If we sum the partial energy-momentum tensors over all constituents of the mixture and compare the resulting 
equation with the energy-momentum tensor of the mixture 



we can identify the quantities of the mixture 

P (a « = Ep^> P = EP„, «7 = 5> ot ne = ^n a e a , 



(3.10) 

(3.11) 
(3.12) 



Note that according to (3.12b) the heat flux of the mixture is a sum of the partial heat fluxes and a term which 
represents the transport of heat due to diffusion. 



IV. GRAD'S DISTRIBUTION FUNCTION 



By supposing that all constituents are at the same temperature T - which represents the mixture's temperature - 
we may describe a relativistic gas mixture by 13r + 1 basic scalar fields: 





r scalar fields; 




3 scalar fields; 


J a ' 


3(r — 1) scalar fields: 


T, 


1 scalar field; 


TU a , 


r scalar fields; 




3r scalar fields; 


< Pa , 


5r scalar fields. 



Unlike the non-relativistic theory of gases, the temperature in the relativistic theory cannot be defined in terms 
of the distribution function. It appears as the reciprocal of the integrating factor of the Pfaffian form de — pdn/n 2 , 
which identifies the potential of the Gibbs equation in equilibrium with the entropy per particle s such that ds = 
^ (de — -^dn) . Furthermore, in this work we are interested in process close to equilibrium where J", w a , q" and pi"^ 
are considered as small quantities. 

In terms of the 13r + 1 basic fields the distribution function of constituent a can be represented as the following 
polynomial function of the momentum four-vector: 

/„ = /<°> (1 + A a a Paa + Af Paa PaP ) , (4.2) 

where A" and A" 13 are tensorial coefficients to be determined from the definitions of the partial particle four-flow 
(2.8a) and partial energy-momentum tensor (2.8b). Furthermore, fa°^ is the Maxwell- Jiittner distribution function 



(4.3) 
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Above k is the Boltzmann constant, K 2 (( a ) denotes the modified Bessel function of second kind defined through the 
integral 

'cy_r(i/2) 



K n (C) 



-Cv I 2 -i \ n-1/2 j 



(4.4) 



r(n + l/2) A 

and Ca = m a c 2 /kT is a parameter which represents the ratio of the rest energy m a c 2 of a relativistic a— particle and 
the thermal energy of the mixture kT. If ( a 3> 1 the constituent behaves as a non-relativistic gas, while when Ca <C 1 
as an ultra-relativistic gas. 

In order to determine the tensorial coefficients we decompose the momentum four-vector p" and the tensorial 
coefficients A" and A"^ as 



Pa = &pP'a + \u a {pP a Up) , A a a = \,U a + A?A^, 



1 



Af = \ 2 U a U + A 3 ?7 Q/3 + (A*£/^ + A^C/ Q ) + A^A^Af . 



(4.5) 
(4.6) 



Here we have introduced new coefficients: the scalars Ai, A2 and A3, the four- vectors A" and K 2 and the symmetric 
and traceless second-order tensor A^ 7 ^. 

Hence, the distribution function (4.2) can be written as 

fa = /i 0) (l + \i P a a U a + \ 2 p a a piU a Up + X 3 m 2 a c 2 + A?j£A a/J 

+Klp a a P^ ia U p + A^p^ a A a ,A p5 ) , (4.7) 

thanks to (4.5) and (4.6). 

Now the insertion of (4.7) into the definition of the partial particle four-flow (2.8a) and of the partial energy- 
momentum tensor (2.8b) and integration of the resulting equations leads to the determination of the lambda tensorial 
coefficients and it follows Grad's distribution function for the constituent a, namely, 



1 — 5G a ( a — Ca + GaCa 



Pa 



p a 20G a + 3Ca - UGlCa - 2G a Q + 2GIQ 



^-jU a u p pyp 



3Ca 



•2 , ^3 



6G a + Ca — G 2 a Q a JT a 15G a + 2C, a ~ GG 2 Ca + 5G a Ca + Ca — ^aCa 



l a C 2 1 — 5G a Ca — Ca + ^aCa 



1 — 5GaCa — Ca + ^aCa 



Ca 



Pa Ca + 5G a — G 2 Ca 



G a 



>Pa 



Pa{a(S) Ca 1 



Pa 2Ga m 2 a c 



(4.8) 



Above we have introduced the abbreviation G a = K 3 (( a )/ K 2 (( a ). 



V. CONSTITUTIVE EQUATIONS 



In order to determine the constitutive equation for the third-order moment tensor of constituent a we insert Grad 
distribution function (4.8) into its definition (2.8c) and integrate the resulting equation, yielding 

= (n a C la + C 2a w a ) U a U^ + j (n a m 2 - n a C la - C 2a w a ) (rf^lP + rf^ 

+^U a ) + C 3a ( V a ^Z + , 7 « 7 qf + ^q«) - |c 3a (U a U f3 ql + U a U^ f3 a + C/% 7 q«) 

+G 4 a (p^lP + pi^U? + pi^U a ) + C 5a (rf#r a + rf^f a + rffj^j 

+^ (u a u p r a + u a wf a + u^vsz 



The scalar coefficients Gi a through G6a are given by 



Gla = -^(Ca+6G a ), 
Sa 



G3a 



m a Ca + 6G a ~ GK a 



C 2a 



Ca Ca + 5G a — G^Ca 

6m a 2C a 3 - 5Ca + (19Ca 2 - 30)G a - (2C| - 45Ca)G 2 - 9C a 2 G a 



C 2 Ca 



20G a + 3Ca - 13G2C a - 2C Q 2 G a + 2Ca 2 G a 



G4a = 7=r7-(Ca + 6G a ), G 5a = ~m a kTG a , G 6a = m a kT (£ a + 6G a ) 

(-raCa 



(5.1) 

(5.2) 

(5.3) 
(5-4) 
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Furthermore, the pressure p a , the energy per particle e Q and the enthalpy per particle h a of constituent a read 

p a = n a kT, e a = m a c 2 ( G a - 7-) , h a = m a c 2 G a . (5.5) 



Ca 

For the determination of the production terms (2.11) and (2.12) it is necessary to introduce the total momentum 
four- vector P a and the relative momentum four- vector Q a defined by (see e.g. [7]) 

P«=K+K, P /a = P T+p' b a , Q a = P a a ~P^ Q' a =p'«-p' b a - (5.6) 

From the above equations together with the momentum four- vector conservation law it follows that 

P a = P' a , P a Q a ^{m 2 a -m 2 )c 2 1 Q 2 = P 2 -2(m 2 a +m 2 )c 2 1 (5.7) 

where P 2 = P a P a and Q 2 = —Q a Q a denote the magnitudes of the total and relative momentum four-vectors, 
respectively. The inverse transformations of (5.6) are 

pa y^)Q pa f\a pa D la P a D ,OL 
„« _ _ . r> a - — — r! a - — 4- -1— r>' a - — -1— K «1 

Pa _ 2 + 2 ' Pb ~ 2 2 ' Pa ~ 2 + 2 ' Pb ~ 2 2 ' 1 j 

The relative momentum four-vector can be written as 

p2 + p 



Q a = ° 02 " ^ + -\/ P4 - 2P2 (^a + + K - m2 ) 2c4 ' (5-9) 



where k a is a spacclike unit vector orthogonal to P a , i.e., k a P a = 0. 

From now on we shall restrict ourselves to the case where the rest masses of the particles of the constituents are 
not too disparate so that m b ~ m a (l + e), (V6 = 1, . . . ,r ^ a) with e denoting a small quantity. In this case we have 
that m 2 — m 2 rj (m a + rrib)e and the term (m 2 — m 2 ) 2 in (5.9) can be neglected. In this case the relative momentum 
four- vector can be approximated by 

Q a — (TO °~^ )C2 P" + Qk°. (5.10) 

and the invariant flux (2.4) reduces to 

PQ 
2 



F ba = Vtofca) 2 - m>gc* - (5.11) 



For the determination of the production terms (2.11) and (2.12) we follow [17, 20] and introduce the invariant 
differential elastic cross-section for Maxwellian particles which can be written as 

= £q f{ ^ (5 ' 12) 

where J-{&) is an arbitrary function of the scattering angle O. By taking into account the expression (5.11) and (5.12) 
we obtain that 

FbaCTab = ^(e). (5.13) 

4c 

This methodology is similar to the one employed in the non-relativistic case and the simple expression (5.13), which 
does not depend on the relative momentum four-vector, will permit us to integrate the production terms without the 
use of Grad's distribution function. 

The calculation of the production terms is schematized in the Appendix A. Here we list only the final results 



p " = e {( t « a - T «y a ) n »<* - { T i 3a - T »y a ) N aa \ , (5.i4) 

6=1 
r (- 

r> a P _ ___ K 'ST y ) T«^7 AT r T a ^ AT j_rpa/3rpj rpaprp^ 

~ 2? ' I JVa 7+ J a J b 7 1 b J a 7 



b=l 



5) ^TJ 7 + 4 (r^T^ + If I*) + (ml + m 2 )c 2 (jV?< + N?N? 



^N bl +T^-<N ai )^2r 1 aP 



T7 s T h ^ - ( i-*- + ^ ] TIT* 



(5.15) 
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where we have introduced the following abbreviation for the integral 

B n = I T{Q) (1 - cos" 6) sin 6d9. 
Jo 



(5.16) 



We note that with this methodology we were able to write the production terms (5.14) and (5.15) in terms on the 
moments of the distribution function N", and T"^ 7 , without the knowledge of the distribution function. 



VI. LINEARIZED 13r+l FIELD EQUATIONS 

Once the constitutive equations for T"^ 7 , P" and P"^ are given in terms of the 13r+l fields (4.1) through (5.1), 
(5.14) and (5.15), respectively, we can obtain from the balance equations (2.10) the corresponding 13r+l field equa- 
tions. Here we are interested in the field equations that are linear in the non-equilibrium quantities J", w a , q", pf^ 
and their derivatives. Hence, the linearized field equations for the partial fields particle number density n a , four- 
velocity U a and temperature T read 



Dn a + n a V a C/ a + V Q J^ =0, 



V a (p 



1 



Dq a = 0, 



n Cu Dr+pV Q t/ a + V Q q Q = 0. 



(6.1) 
(6.2) 
(6.3) 



The balance equation of the particle number density (6.1) follows from (2.10a) by using (B2) of Appendix B. Note 
that we have introduced the operators D and V Q defined through the relationships 



D = U a d ai V Q = Aidp 



so that d a = -^rU a D + V Q 

c z 



and 



£/ a V Q = 0. 



(6.4) 



For the mixture the balance equations of the four-velocity (6.2) and of the temperature (6.3) follow from (2.10b), (B3) 
and (B5) of the Appendix B and by summing the resulting equation over all constituents. The former refers to the 
projection of this equation by A a ^, while the later is the projection with respect to U a . In (6.3) it was introduced 
the heat capacity per particle at constant volume of the mixture 



(6.5) 



We get the balance equations for the partial diffusion fluxes J" from (2.10b), (B3) and (B5) of the Appendix B by 
taking the projection A« of the resulting equation and by subtracting the rth equation from the ath equation. This 
is necessary in order to obtain 3(r — 1) independent scalar balance equations. Hence it follows 



D 



n Q h a 



VppW - V" (p a + w a ) + -^Dq«l - -f- [v^> + iDq= 
^—By V ( [h Q + h & - 3kT] (n 6 J£ - n J?) + n b q« - n a q«) 

^—B x V ( [h 6 + h r - 3fcT] (n r J£ - n b J?) + n r q£ - n b q«). 
n r h r ^ [ J 



-V a (p r + m r ) 



(6.6) 



The balance equations for the partial non-equilibrium pressures zu a , heat fluxes q" and pressure deviator tensors 
p^"^ are obtained from (2.10c), together with (B4) and (B6) of the Appendix B as follows: first the projection U a Up 
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lead to the balance equations for the partial non-equilibrium pressures w a : 



-Dw n 



C*5a + Cg a 



Dn a - ^n a C' la DT - 5-^V 7 q^ + f(rn 2 a + 5<7 lo )V 7 EP 
r 

V 7 J a = "TT^ 1 X! [ nfc (^^a + 6e b ) W a - (l a (c 4 C 2b + 6e a ) G7 b ] 



6=1 



6=1 ^ 



C 2a - 3 5 1 e 6 + 6 (e 6 - 3fcT) —2- + 



m a m b 



+6fcT 
+6fcT 



n a tu b 
3n a n b 



c 4 C 2b - 3 5 1 e a + 6 (e a - 3fcT) —2- + 



2m& 2m a 



2mj 2m Q 



(mj - m\) (e a e b - 3(fcT) 2 ) - 3fcT (m 2 e a - m 2 e b ) 



-kT (m 2 a e b - m b e a ) 



9 9 f ' 
™a m b 



(6.7) 



where we have introduced C( a = dCi a (( a )/d( a . Next, the balance equations for the partial heat fluxes q" follow from 
the projection A"C/^: 



5C 3a Dq^ 



6 



[ml - Ci a ) V Q n a + ^n Q C( Q V Q r - C 2a W a m a 



2 ' 

(C 5a + C 6a ) DJ^ - - (m 2 n Q + 5n Q C la ) DU a = ~Bi £ \ c 2 n b J« 



G 



6=1 



6 

hh , s (m 2 , — ml)c 2 

+4(e a - 3fcT) + V a fe; 



c 2 nj b 



+ Cgb 



C-5a + C*6a 

(™ 2 - C la )c 2 



6 



+ n 5 (5c 2 C 3a - e & + 3fcT)q« - n a (5c 2 C 3 6 - e a 



6=1 



10C 3a c 2 + 4(e b - *T) - 5 1 (e 6 - 3fcT) 



+n a 



10C 3b c z + 4(e a - fcT) - ( 5 1) (e a - 3fcT) 



77! 



q? + 2c^n b 



C5a + C*6a 



(m 2 -C lb )c 2 A.m 2 ^ e b - 3fcT h a 



2^(e b -fcT) 



(m 2 +m 2 )c 2 



-* n 



+2c 2 n a 



(m 2 a + m 2 )c 21 



(m 2 a -C la )c 2 
6 



ml \ e Q - 3/feT h b „ h 



2^(e a -fcT) 



(6.8) 



Finally, the projection A" Ajj — A 7< 5A"^/3 imply the balance equations for the partial pressure deviator tensors pi' 



<«/3>. 



C 4a Dp<^ + 2C 3a V< Q q^ + 2C 5a V< Q J^ + yn a (m 2 - C la ) V< a £/« 

= ~Bi E |n6(c 2 C 4a + e b - 3fcT)p<^> - n a (c 2 C 4b + e a - 3fcT) p< a/s >| 

C 6=1 ^ J 

6=1 



2c 2 C 4a + 5 1 (e 6 - 3fcT) + 8fcT 

mt 



p <-0 



+n a 



m| 



2c 2 C 4b + 5 1 (e„ - 3kT) + 8kT 



(6.9) 



Hence, the system of equations (6.1), (6.2), (6.3), (6.6), (6.7), (6.8) and (6.9) compose 13r+l linearized field 
equations for the fields (4.1). 
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VII. SIX-FIELD THEORY FOR A BINARY MIXTURE 

In this section we shall restrict ourselves to a binary mixture characterized by the six scalar fields of particle number 
densities rq, n 2 , four-velocity U a and temperature T, whose balance equations (6.1) - (6.3) are written as 

Dm + mValT" + V a J? = 0, Dn 2 + n 2 V Q C/ Q + V Q J^ = 0, (7.1) 

^DU a +V, 3 p {a0) - V° (p + zu) + ^Dq" = 0, ncDT + pV l7 a + V a q a = 0. (7.2) 

In this case the pressure deviator tensors pj"^\ p 2 "^\ the non-equilibrium pressures w\, nj 2 , the heat fluxes q" , q 2 
and the diffusion flux J" = — J 2 are no longer variables, but constitutive quantities. To determine these constitutive 
quantities we shall rely on the remaining 21 scalar equations (6.6) - (6.9) and a method akin to the Maxwcllian 
iteration procedure, which is often used in kinetic theory of gases. In this method the equilibrium values of the 
constitutive quantities - namely, p x = p 2 a ^~0, w\ = zu 2 — 0, q" = q 2 = and J" = — J 2 = - are inserted on 
the left-hand sides of the remaining 21 scalar equations and the first iterated values are obtained from the production 
terms, i.e., from the right-hand sides of these equations. 

Following the methodology described in the Appendix C we get from: 

(i) two equations that follow from (6.9) 

p< Q/3) =2 W V<°1^\ P { 2 P) =2^ 2 V< q l^; (7.3) 

(ii) two equations that follow from (6.7) 

w x = -mV a U a , m 2 = -Ti 2 V a U a ; (7.4) 

(iii) one equation and two equations that follow from (6.6) and (6.8), respectively, 

^ = -^=V 12 d%+V T V a T, q? = AiVT+X^d?, q 2 = A 2 V a T + 2? 2 d". (7.5) 

Equations (7.3) and (7.4) represent the constitutive equations of a relativistic Newtonian fluid - also known as the 
Navier-Stokes law - and the associated transport coefficients fi\, \i2 and 771 , r\2 are identified as the coefficients of shear 
and bulk viscosities, respectively. In (7.5) we have introduced the generalized diffusion forces 

d " = h ( yQpa_ ^r vap ) ' a ^' 2 ' (7 - 6) 

that are restricted by the constraint Y^ a =i = 0> so that d 2 = — d", and the relativistic temperature gradient 

V Q T = (V*T - ^V"p^) , (7.7) 

which in the non-relativistic limiting case reduces to the temperature gradient V Q T. Hence, we may identify (7.5a) 
as the generalized Fick's law and (7.5b,c) as the generalized Fourier's law. The corresponding transport coefficients in 
these equations are recognized as the coefficients of diffusion 2?i 2 , thermal-diffusion T>t, thermal conductivity Ai, A 2 
and diffusion-thermal T>\ , 2? 2 . 

Now we are ready to obtain the constitutive equations necessary to convert the system of equations (7.1) and (7.2) 
into a system of field equations for the six fields of partial particle number densities m, n 2 , four- velocity U a and 
temperature T. Indeed, by using the definitions of the pressure deviator tensor, non-equilibrium pressure and heat 
flux of the mixture given by (3.11) and (3.12) we get 

p^ = 2fiW {a U p) , w = -riV a U a , q Q = A'V Q T + 23d?. (7.8) 

Here the transport coefficients of shear viscosity fj,, bulk viscosity 77 and thermal conductivity A' of the mixture and 
the diffusion-thermal coefficient T> read 

H = Hi+H2, i] = t] 1 + t]2 7 A' = Ai + A 2 + (hi — h 2 ) T>t, (7.9) 
V = V 1 +V 2 + (h 1 -h2)V 12 . (7.10) 

We call attention to the fact that the true thermal conductivity of the mixture is defined as the ratio of the heat flux 
and the temperature gradient when there is no diffusion , i.e., when J" = — J 2 = 0. In this case the we get from (7.5) 
that 

df = -|^V Q T, hence q Q = AV"T, (7.11) 
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where the true thermal conductivity of the mixture is given by 



T>T 

A = Ai + A 2 — — -- (T>i + V 2 ) . 

V\2 



(7.12) 



The constitutive relation for the diffusion flux (7.5a) together with the ones for the pressure deviator tensor, non- 
equilibrium pressure and heat flux (7.8) imply into the desired field equations for the six fields (ni, ri2, U a ,T), when 
they are inserted into the system of equations (7.1) and (7.2). 

The expressions for the transport coefficients, even in the order up to e, are too large to be given here. Bellow we 
present only their expressions in the non- and ultra-relativistic limiting cases. 

In the non-relativistic limiting case the thermal energy of the gas kT is smaller than the particle rest energy mic 2 , 
so that d> 1 and we obtain 



£>12 = 



kT 



(gi + ff 2 )n + 2n 1 ff 1 4gi + B 2 
2n(2i3 1 +S 2 ) 6 2(26! + B 2 ) 



1 



B 2 (26m + 6n 2 ) + Bf (m + 3n 2 ) + BiB 2 (13m + lln 2 



V 
f 



2n(2B 1 + B 2 )(4B 1 +B 2 ) 
5(Bi+B 2 )fcmn 2 fl 5B 2 + Bi£ 2 -B 2 1 



47rmm 2 B 2 (2Bi + B 2 ) {(i (Bj + B 2 ){2B 1 + B 2 ) Q 2 



kTh\ 



7TTOin B\B 2 



2kT 



1 



1 



3ttB 2 
5kT 1 

5k 2 T 



Bi+B 2 - 
1 - — 



Bi(6Bi + 7B 2 ) 1 



4B! + 2B 2 Ci 
n 



l-2^ e - 25 



n 



2Ci 



1-3^6 



2— e 



1 



■ o 



(l 2 



Ci 3 



'd 2 

1 



A = A' = 



27rmi B 2 



n, 39 



l-3^e 



(7.13) 
(7.14) 
(7.15) 
(7.16) 
(7.17) 
(7.18) 



In the ultra-relativistic limiting case the thermal energy of the gas kT is larger than the particle rest energy mic 2 , so 
that and we have 



V l2 = 

V T 
V 

/' 
V 



6c 2 (2B 1 +B 2 ) 
7rnBi(26Bi + 11B 2 ) 



1 - 



2(llBi + 5B 2 ) 
3(26Bi + ilBal 



1 



4B 2 (15m + 7n 2 ) + BiB 2 (55m + 29n 2 ) + B|(13m + 7n 2 ) 
2n(llBi + 5B 2 )(2Bi+B 2 ) 
c 2 (6Bi + llB 2 )nin 2 
_ 37rn 2 TB 2 (26Bi + 11B 2 ) 

kTc 2 26B 2 + hhB x B 2 + 26B 2 2 



1 



7rnB 2 

AkT 
5ttB 2 

kTn 1 
1087rnB 2 

4c 2 k 



78Bi + 33B 2 



80 



1 



, n 2 



eC 1 2 +0(e 2 ,Ci), 

Cb + 0(e 2 ,Ct), 

Ci 



(7.19) 
(7.20) 



e + 7 + In 



1 



<D(e 2 Xt) 



CI 1 



10n 2 (22Bi +29B 2 ) 



A = A' = 



3ttB 2 



1 - 



n(2Bi 
1 



B 2 



n 2 
2— e 



e }+0( e 2 ,d 5 ), 

o( £ 2 ,Ci 3 ), 



(7.21) 
(7.22) 
(7.23) 



From the expressions of the transport coefficients (7.13) - (7.23) we note that: 



1. they are valid up to the e-order, which was the approximation used when the masses of the particles of the 
constituents are not too disparate, so that we have written m 2 = mi(l + e) with e being a small quantity; 

2. in all expressions relativistic corrections are given in terms of the parameter d; 

3. the coefficients of thermal-diffusion T>t and diffusion-thermal T> are of order e; 
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4. the thermal-diffusion coefficient T>t given by (7.14) is of relativistic order, i.e., it vanishes in the non-relativistic 
limiting case. This is well known result in the literature that the thermal-diffusion coefficient vanishes for a 
mixture of Maxwcllian particles (see e.g. [25, 26]); 

5. the diffusion-thermal coefficient T> is also of relativistic order, since its expression (7.15) depends on the rest 
energy = ra\c 2 ] 

6. within the e— order the coefficients A and A', associated with the thermal conductivity of the mixture, coincide; 

7. for a single component - i.e., when X2 = - the transport coefficients of shear and bulk viscosities and thermal 
conductivity in the non-relativistic and ultra-relativistic limiting cases reduce to the ones given in [17, 20]. 

In the case of the particles have the same rest masses m\ = m 2 the coefficients of thermal-diffusion T>t and 
diffusion-thermal V vanish and the diffusion coefficient becomes the self-diffusion coefficient T>u whose expression is 
given by 

= -2^(25. + Ba )(10G 1+ Ci-G;Ci) > where (7 . 24) 

A = rrnB 1 {2B 1 [G?C 2 (GiCi - 8) + Gi(30 + 8(f) + £(3 + C?) + 2G?Ci(l - (!)] 

+B 2 (30G 1 +3Ci - 13G?Ci)}- (7-25) 

When mi = m 2 the character of the mixture is exclusively owing to the difference of the particle number densities 
of the components. Furthermore the coefficients of thermal conductivity, shear and bulk viscosities of the mixture 
reduce to 

= 2fcc 2 C(C + 5G - G 2 Q 2 = 2fcTG 2 C 

7r£ 2 (C + 10G-CG 2 ) ' ^ ^? 2 (2G + C + 2G 2 C)' 1 ' ' 

_ kT((2QG + 3C - 13G 2 C - 2GC 2 + 2G 3 C 2 ) 2 
V ~ 3tt6 2 (1 - 5GC - C 2 + G 2 C 2 ) 2 (C + 5G - CG 2 ) ' ( ' 

which are the same as those of a single gas [17, 20]. The constitutive equations in the case where the masses of the 
particles arc identical read 

J? = -Jf = £>n df , p <Q,3> = 2 A tV (Q [/' 3) , w = -vV a U a , q" = AV Q T. (7.28) 

In the next section we shall investigate the solutions of the six-field theory concerning the propagation of forced 
and free waves in a relativistic binary mixture and for that end we shall need the values of the transport coefficients: 
diffusion T>i 2 thermal-diffusion T>t, thermal conductivity A', diffusion-thermal T>, shear viscosity \x and bulk viscosity 

V- 

In Figures 1-3 arc plotted the dimensionlcss coefficients 

Vt 2 =V 12 *f, V* T =V T ^, A*=A'0, (7.29) 

^* „ $2n „ B 2 * B 2 , , 

V =V Wc^ " =/i fcT' 11 =7] kT' (7 - 30) 

as functions of the parameter £ = m\c r/kT. In these figures it was consider that the concentration of the constituent 
labeled with the index 1 was 60%, i.e., xi = ni/n = 0.6 with x 2 = 1 — xj., and that the mass of the constituent 2 
was 10% larger than that of the constituent 1, so that e = 0.1. Furthermore, the ratio of the integrals B 2 jB\ was 
taken equal to the ratio of the integrals that appear in the theory of non-relativistic Maxwellian particles, namely, 
A 2 (5)/Ai(5) (see, e.g. [25, 26]). Hence, we have adopted the value B 2 /B x = 0.436/0.422. 
We infer from these figures for the dimensionlcss coefficients that: 

1. the values of the diffusion 2?* 2 , shear viscosity /x* and thermal conductivity A* in the non-relativistic limiting 
case are smaller than those in the ultra-relativistic one; 

2. the coefficients of thermal-diffusion T>^ and bulk viscosity 77* vanish for very small and very large values of £1; 



3. the diffusion-thermal coefficient T>* is very small in the ultra-relativistic limit and tends to a constant value in 
the non-relativistic limit. 



13 






FIG. 2. Dimensionless coefficients as functions of fi: thermal conductivity A* (fig. 2a) and diffusion-thermal T>* (fig. 2b). 



VIII. FORCED AND FREE WAVES IN A RELATIVISTIC BINARY MIXTURE 

As was pointed out in the last section,, the insertion of the constitutive relations (7.5a) and (7.8) into the balance 
equations (7.1) and (7.2) leads to system of field equations for the determination of the six fields of partial particle 
number densities n a (a = 1,2), mixture four- velocity U a and mixture temperature T. Let us then look for solutions 
of form 



U' 



n a = x a n + n a exp (iqx — iuit) , T = To + T exp (iqx — iujt) 
c Uexp (iqx — not) 



yl — [{v/c) exp (iqx — ttot)] 2 \J 1 — [(v/c) exp (iqx — not)] 



= ,0,0 , 



(8.1) 
(8.2) 
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1.5e-3 



le-3- 



5e-4- 




FIG. 3. Dimensionless coefficients as functions of £1: shear viscosity /i* (fig. 3a) and bulk viscosity rf (fig. 3b). 



which represent small perturbations about an equilibrium state characterized by constant partial number densities 
x a no, constant temperature To and vanishing value of the velocity in the longitudinal direction parallel to the x axis. 
Besides, the amplitudes of the perturbed fields - i.e., the overlined quantities - are small so that only linear deviations 
from equilibrium need to be take into account, while oj and q are angular frequency and the wavenumber of the 
perturbations, respectively. 

In the linear case, we can easily verify that the operators D and V" become D = d/dt and V Q = (0, —d/dx l ). Hence, 
the insertion of expressions (8.1) and (8.2) into the field equations leads to following system of algebraic equations for 
the amplitude of the perturbations 



/Mn M u M13 M w 

M 2 i M 22 M 23 M 24 

M 3 i M 32 M 33 M 34 

V M 4 i M 42 M 43 M 44 



( ni/no ^ 






n 2 /n 







v/c 







{ T/T J 




U 



where the elements of the matrix are given by 

Mn = + 1 

Mi 3 

M 2 i 



x 2 (1 + eaiXi)X>* 2 - — ^ 
n* 



xi(l-eaix 2 )Hia + T z 



Xig*, M14 = i 
x 2 (1 + eaiXi)X>i2 



eaixix 2 l?i 2 + — (h* - 1) 

V, " 



M 22 = a;* + i 
M 24 = -M u , 



xi(l-eaix 2 )2?* 2 + ^ 



M 23 = -x 2 g* 



M 3i 



-g* + 1 



x 2 (1 + eaiXi)T>* 



M 32 
M 33 



"9* 



— +Xi (1 - eaix^rJ' 



g*w*, 



(8.3) 



8.4) 
8.5) 
8.6) 
8.7) 
8.8) 
8.9) 
(8.10) 
(8.11) 
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M 34 
M A2 



= -q*-i 



(h* - 1) + eaiXix 2 X>* 



x 2 {l + ea 1 x 1 )V* 



A* 

xi (1 - eaix 2 ) V* + — 
n* 



q*uj* 



M 43 = 



Maa = C* + 2 



eaixix 2 2?* 



(8.12) 
(8.13) 
(8.14) 
(8.15) 



In the above equations we have introduced a mean free time r and starred quantities which correspond to dimensionlcss 
coefficients. They are given by 

1 



<1* 



qcr, 



= LOT, 



h* = 777; = C1G1 (1 + e CK1X2) , 



gfCi - Ci - 4gi 
Gi 



(8.16) 
(8.17) 



< = J = Ci 2 + 5CiGi - Ci 2 G? - 1 + ex 2 Ci [2d 2 G! (l - Gj) + Ci (13G? - 3) - 20G X ] 



(8.18) 



The system of algebraic equation (8.3) has a non-trivial solution if the determinant of the matrix of the coefficient 
vanishes. This condition leads to a dispersion relation which can be used to study the propagation of sound waves, as 
well as the dynamical behaviour of small disturbances induced by spontaneous internal fluctuations. For sound waves, 
the dispersion relation is solved by taking the angular frequency as a real input variable. As a result a wavenumber is 
found which in general is complex. In this case, the phase velocity v p and the attenuation coefficient a are defined by 



and 



5.19) 



On the other hand, for the eigenmodes, we consider the wavenumber as real input and solve the dispersion relation 
to obtain the complex angular frequency. The real part of the angular frequency 5ft(w) gives the oscillation frequency 
of a small internal perturbation with wavenumber q, while its imaginary part $s(u)) describes the decay in time of the 
oscillation amplitude. 

A. Acoustic Solution in the Low-Frequency Limit 

In the low-frequency limit, the acoustic solution of the dispersion relation can be determined by expanding the 
reduced wavenumber q* in power series of the reduced oscillation frequency cj* as 



<7* = ciq to* + a\ cj„ + 



(8.20) 



where 00, 01, • • ■ are complex expansion coefficients. By inserting the expansion (8.20) into the dispersion relation 
and equating equal powers of the reduced oscillation frequency w* we get 



ao = ±1 



5.21) 



dl = I < — U + 77 

2j\r 1 



A* +7 [(7-l)a^-l] 2?*a lXl x 2 e 



5.22) 



where 7 = (c* + l)/c* is the specific heat ratio. 

The non-relativistic limit for the phase velocity and the attenuation coefficient in the low-frequency regime can be 
derived from the above expressions when we assume that the thermal energy of the gas mixture is smaller than the 
particle rest energy, i.e., when C,\ 1. Hence, we have 



J = ±jl--^(l-x 2 e) + .. 

av o 7 , . f, 15 ,„ . 
u) 157r I 28^i 



5.23) 
5.24) 
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where vq denotes the adiabatic sound speed for a non-relativistic binary gas mixture: 

v a = ./|^£(l_ X2 e) (8.25) 
V 3 mi 

At this point, it is important to mention that in the derivation of the above expressions for the phase velocity and 
the attenuation coefficient we have assumed that the masses of the particles of the constituents are not too disparate 
so that e< 1. 

When Ci — > oo, we verify from expressions (8.23) and (8.24) that the phase velocity reduces to the adiabatic sound 
speed for a non-relativistic binary gas mixture, while the reduced attenuation coefficient avo/tu varies linearly with 
the reduced frequency. 

On the other hand, the ultra-relativistic limit for the phase velocity and the attenuation coefficient follows when 
(i<l, i.e., when the thermal energy of the gas mixture is larger than the particle rest energy. In this case, we have 

v v If (1 + 2x 2 e) >9 1 , 

^=±_|l-L_±! C? + ...j ( 8 . 26 ) 

££=±^ M /i + £>(£)}. (8 . 27) 

When d — > we note that the phase velocity is of order of the speed of light v p = ±c/ \/3, while the reduced attenuation 
coefficient, as in the non-relativistic limiting case, also varies linearly with the reduced oscillation frequency. 

If we consider x 2 = in (8.23) - (8.27) we get the phase velocities and the attenuation coefficients for a single 
relativistic gas with Maxwellian particles. The expressions for the phase velocities (8.23) and (8.26) when x 2 = are 
the same as those in the work [27]. However, the attenuation coefficients (8.24) and (8.27) agree with those in [27] 
only in the first term without the corrections in the parameter £i. This fact can be easily understood, since here 
we are considering Maxwellian particles while in [27] a kind of hard-sphere potential was used. Due to the fact that 
the attenuation coefficient depends on the transport coefficients, which are functions on the type of potential, the 
corrections in the parameter £i should not agree. 



B. Eigenmodes in the Small Wavenumber Limit 



By solving the dispersion relation to obtain the reduced oscillation frequency w* as a function of the reduced 
wavenumber we verify the existence of four longitudinal eigenmodes. Two of these hydrodynamic eigenmodes are 
sound modes describing sound propagation in opposite directions parallel to the wavenumber, while the other two 
eigenmodes are non-propagating modes that gives rise to purely diffusive effects. Expressions for these eigenmodes 
can be derived in the small wavenumber limit by expanding in power series of q# as 



do + a\ q* 



a 2 <?* 



(8.28) 



By inserting expansion (8.28) into the dispersion relation and equating equal powers of q* we obtain for the sound 
modes: 



a = 0, 



ai = ±\ — 



1_ 



a 2 = - 

and for the diffusive modes 
a = 0, 

a2 = 



2 7 [3 



(7 -I)] 5 



K-(t-i)] 



V*a.i xix 2 e 




A* ± 



T>* 



(7-1) 



7 



(7- 1) 
•4— -V*V* T 

1 



1/2" 



(8.29) 
(8.30) 

(8.31) 
(8.32) 



As mentioned in Section VII the thermal-diffusion coefficient T>t and the diffusion-thermal coefficient T> are of order 
e so that the products V*e in (8.30) and V*V^p in (8.32) can be neglected. In this case, the two sound modes can be 
written as 

[7-M7-1)] 2 



LOT = ±. 



l iqCT) - 1 ^ 



-J 1 



V 



h*7 



A* {qcrY + 



5.33) 
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and the two diffusive modes as 



(7- 1) ^2 , 

LOT = —I A {qcT) + 

7 



and 



LOT = —iD* 12 {<1Ct) 2 + 



5.34) 



We may observe from (8.34) a decoupling between entropy fluctuations and concentration fluctuations. A similar 
result is observed in a non-relavistic binary gas mixture of Maxwellian particles [28] , where thermal-diffusion is zero 
and the non-propagating eigenmodes are related to thermal diffusivity and mass diffusion processes, respectively. 

In the non-rclativistic limit, the two sound modes read 



-^(l-x 2 e) 



7 



33 



(W)-*— U- — (l-x 2 e) + ... \(qv Q TY + .. 



while the two diffusive modes become 



LOT 



5tt 

57T 



(1 



2Ci 
1(3 



- x 2 e 
B) 



15tt 



(qv T) 2 
(l-2x 2 )e 



7( 



1 (4 + B) 
'2(2 + 6) 



1 



2(2 + 6) 
26 - 36x 2 + (13 - 14x 2 )£ 



B 2 



2(2 + B){A + B) 



1 



+ ..A (qv T) 2 + 



5.35) 



(8.36) 



(8.37) 



Above we have introduced the abbreviation B = B2/B1 . 

The two sound modes in the ultra-relativistic limiting case are given by 



4 



1 



— 1- 



157T 



(1 + 2X26) 2 

24 Cl + 

(l + 2x 2 e) A2 
5 M + • 



(qcr) 
(qcT) 2 



and the two diffusive modes read 



LOT 



3tt 
6S 



1 



24 
(2 + B) 



(1 



7T (26+116) 



2x 2 e)Ci 2 + 
1 2 (11 



(qCTf 



5B) 



3 (26 + US) 



d 2 + ... HqcT) 



(60 - 32x 2 + (55 - 26x 2 )B + (13 - 6x 2 )6 2 
+ 2(2 + 6)(26+116) 

Note that fluctuations associated with the sound modes propagate with velocities v = $t(uj)/q given by 



±vq<1 



7 

4Ci 



(l-x 2 e) + 



± 



24 



(l + 2x 2 e)Ci 2 + 



(8.38) 



5.39) 



(8.40) 



(8.41) 



in the non-relativistic and in the ultra-relativistic limiting cases, respectively. 

The expression (8.33) for the sound modes and the one for the diffusive mode (8.34a) are the same as those obtained 
in the work [29] for a single relativistic gas. If we set X2 = in (8.33) and in (8.34a) we get the sound and diffusive 
modes of a single relativistic gas with Maxwellian particles. As was pointed out in the last subsection the relativistic 
corrections in the parameter (j should not agree with those in [29], since in the latter work a kind of hard-sphere 
potential was used. 

The frequency and wavenumber regions of validity of the acoustic and eigenmodes solutions derived from the six 
field theory for a binary mixture can be determined by imposing, respectively, the conditions (see [29]): 



and 



qpc 
P 



< 1, 



5.42) 



where fi is the shear viscosity of the mixture in the non-rclativistic limiting case. The first condition tell us that an 
acoustic solution based on the classical hydrodynamic description is valid as long as the molecular collision frequency 
is larger than the wave frequency, while the second one tell us that the eigenmode solution is valid as long as the mean 
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free path of the molecules is smaller than the wavelength of the internal fluctuations. In terms of the dimcnsionlcss 
frequency u>* = ujt and the dimensionless wavenumber = qcr we can rewrite the above conditions as 

Sir Sir 
uj*> —, and g* < — . (8.43) 

Finally, it is important to mention that another possible way to test the range of validity of an extended hydrodynamic 
description for the acoustic problem in rclativistic gases is to consider how the maximum speed of propagation 
approaches the light speed with the increasing number of the moments. For the rclativistic single gas case, a complete 
analysis based on this method was given in the work [30] . 



IX. CONCLUSIONS 



To sum up a mixture of r constituents was analyzed within the framework of Boltzmann equation by using Grad's 
moment method. The mixture was described by 13r + 1 basic fields of four-velocity, temperature of the mixture, 
particle number densities, diffusion fluxes, non-equilibrium pressures, heat fluxes and pressure deviator tensors. The 
13r + 1 field equations were obtained from Grad's distribution function applied to a mixture of Maxwellian particles 
where the rest masses of the constituents are not too disparate. From the system of 13?* + 1 field equations it was 
determined the constitutive equations for a binary mixture with six scalar fields of particle number densities, four- 
velocity and temperature. The generalized laws of Fick, Fourier and Navier-Stokes were obtained and the transport 
coefficients of diffusion, thermal-diffusion, diffusion-thermal, thermal conductivity and bulk and shear viscosities were 
determined. Explicit expressions for these coefficients were given in the non-relativistic and ultra-relativistic limiting 
cases as well as their graphs in terms of a parameter which represents the ratio of the rest energy of a particle and the 
thermal energy of the gas mixture. An analysis of forced and free waves for the binary mixture was also performed. 
The phase velocity and attenuation coefficient were determined for the acoustic solution in the low frequency limit. 
It was also shown that there exist two sound modes describing sound propagation in opposite directions parallel to 
the wavenumber, and two non-propagating modes related to purely diffusive effects. 
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Appendix A: Evaluation of the production terms 



Let us evaluate the production term Pf. By using the relationships (5.8) and the representation (5.13), the 
production term (2.11) becomes 



6=1 •* 
r „ 

= o£ / Q( - k " 3 - kP )fahF b a<Ja b dCt 
6=1 •* 



d 3 p b d 3 p a 

PbO PaO 

d 3 p b d 3 p a 

PbO PaO 



(Al) 



If we write the element of solid angle as dfl = sin OdOdQ where G and are the spherical angles of k'^ with respect 
to VP and use the following result 



p27T fK p 7T 

/ / (k' f3 - k fj ) sin SdSd® = -2ir / 1^(1 - cos 9) sin 9d9 
Jo Jo Jo 

the production term (Al) reduces to 

= E f fQ^fafbF ba * ab dCl d3pbd3pa 
b=l J ° J 



(A2) 
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p2 
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PbO PaO 
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1 — cos 9) sin 9<i9, 



(A3) 
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thanks to (5.10). Now if we take into account (5.13) and transform the variables (P a ,Q a ) into (Pa,Pb) by using the 

d 3 p b d 3 p a 

PbO PaO 



relationships (5.6), we may perform the integrations of (Al) with respect to 4^71777 an d get 



N„ 



C 6=1 

In the above expression B\ is the following integral when n = 1: 

B n = [ J"(9) (1 - cos" 6) sin 9d9. 
Jo 

The same methodology may be used to determine the production term P" . First we obtain from (2.12) 

d 3 p b d 3 p a 
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Next by using the relationship 
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and (A2) it follows that the production term (A6) reduces to 
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By performing the integration with respect to leads to 
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It is worth to call attention that (A2) and (A7) result from the integration in the spherical angles, when the relative 
momentum four-vector after collision is expressed in terms of the one before collision. These are general expressions 
which are valid for all representations of the distribution function, like those that appear in the Chapman-Enskog and 
Grad methods. 



Appendix B: Linearized balance equations 

By introducing the operators D and V a defined by the relationship 

1 



D = U a d a , 



A Q/3 9«, so that d a = -^U a D + V Q , and U a V a = 0, 



(Bl) 
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the linearized left-hand sides of (2.10) read 

d a N« = Dn a + n Q V Q C/ Q + V Q J^, 



DU a + V^p< a « - V" (p a + w a ) + D (q« + u a r a ) 
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u a 
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thanks to the representations (3.2), (3.6) and (5.1). Above we have introduced the notation C' la = dCx a {C, a )/dC, a . 

Again by using (3.2), (3.6) and (5.1) we obtain from (A4) and (A9) the linearized production terms as functions of 
the 13r+l moments of the distribution function: 
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Appendix C: Maxwellian iteration procedure 



In the Maxwellian iteration procedure the equilibrium values of the constitutive quantities pi = p 2 — 0, vd\ = 
w 2 = 0, q" = q 2 = and J" = — i 2 = are inserted on the left-hand sides of the equations (6.6) - (6.9) and the first 
iterated values are obtained from the production terms, i.e., from the right-hand side of these equations. 

First from (6.6) and (6.8) we obtain - by following the above mentioned methodology - the following algebraic 
system of equations for the determination of the diffusion J" and heat q Q , q 2 fluxes: 
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The solution of the above system of algebraic equations leads to the constitutive equations for J", q" and q% : 

j? = -j« = Vvx d? + p T v Q r, q? = AxV^r + Pid?, q j = A 2 v Q r+p 2 d5 i , 

which are linear functions of the relativistic temperature gradient 



V Q T = V Q T- -^V a p, 
nh 



and of the generalized diffusion force 

da 
1 



fcT I Pl nh 



with d£ 



Equation (C4a) represents generalized Fick's law, while (C4b,c) generalized Fourier's law. Here Ax, A 2 , 2?x2, T>t 
T>2 are scalar coefficients. 
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The pressure deviator tensors p[ a ^\ P 2 "^ are obtained from the algebraic system 
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that follow from (6.9) by using the Maxwellian iteration procedure. From the above system of algebraic equations 
we obtain the constitutive equations for the pressure deviator tensors as linear functions of the traceless part of the 
velocity gradient V^ a U^\ namely. 



p< a/J> = 2 m V< Q [/ /3 \ p< a/3> = 2/i 2 V< Q [A 



(C9) 



where /ii,/x 2 are scalar coefficients. The traceless part of the gradient of velocity is defined in the same way as the 
tensor T Q/3 in (3.9). 

Finally, the Maxwellian iteration procedure when applied to (6.7) leads to the following system of algebraic equations 
for the non-equilibrium pressures -071,1172: 
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We note here that the above underlined terms are of order e 2 which must be neglected due to the hypothesis that the 
masses of the particles of the constituents are not too disparate, i.e., m 2 = toi(1 + e) with e being a small quantity. 
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Hence we obtain from (CIO) and (Cll) the linear relationships between the non-equilibrium pressures and the velocity 
divergent \7 a U a : 



ZD! = -i]!\7 a U a , w 2 = -rj 2 V a U a , 



(C12) 



where 771,772 are scalar coefficients. 

Equations (C9) and (C12) represent the constitutive equations of a relativistic Newtonian fluid, also known as 
Navicr- Stokes law. 



Appendix D: Generalized diffusion force and relativistic temperature gradient in the e-order 



In the e-order the generalized diffusion force and the relativistic temperature gradient read 
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